AfRne model of stress stiffening in semiflexible filament networks 



OO 

o 
o 

bjo: 
< 

On 



J.R. Blundell and E.M. Terentjev 
Cavendish Laboratory, University of Cambridge, Madingley Road, Cambridge, CBS OHE, U.K. 

(Dated: August 29, 2008) 

We present a revised theoretical study of the affine assumption applied to semiflexible networks. 
Drawing on simple models of semiflexible worm-like chains we derive an expression for the probability 
distribution of crosslink separations valid at all separations. This accounts for both entropic and 
mechanical filament stretching. From this we obtain the free energy density of such networks 
explicitly as a function of applied strain. We are therefore able to calculate the elastic moduli of 
such networks for any imposed strain or stress. We find that accounting for the distribution of 
cross-link separations destroys the simple scaling of modulus with stress that is well known in single 
chains, and that such scaling is sensitive to the mechanical stretch modulus of individual filaments. 
We compare this model to three experimental data sets, for networks of different types of filaments, 
and find that a properly treated affine model can successfully account for the data. We find that 
for networks of stiffer filaments, such as F-actin, to fit data we require a much smaller effective 
persistence length than usually assumed to be characteristic of this filament type. We propose that 
such an effectively reduced rigidity of filaments might be a consequence of network formation. 

PACS numbers: 82.35.p, 78.20.Ek, 87.19.R 
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I. INTRODUCTION 

Networks of crosslinked or branched semiflexible fila- 
ments are found in many biologically relevant systems. 
These include the networks of filaments such as F-actin, 
intermediate filaments and microtubules that make up 
the cytoskeleton of cells, as well as networks of collagen 
and fibrin found in the extracellular matrix. Looking 
broader, similar semiflexible/rigid networks are formed 
by segments of DNA [1], by aggregated amyloid fibrils 
[2j, self-assembled peptide nanotubes [3] and by carbon 
nanotubes, on their own or dispersed in polymer matrix 
0]. Semiflexible networks of cytoskeleton have interest- 
ing and unusual elastic properties that are thought to be 
crucial to the way in which cells move, function and re- 
spond to their surroundings @, @, 0]- Iii particular, such 
networks tend to be stiffer than conventional polymer 
networks such as rubber, and show a dramatic stiffening 
over modest strains that is absent in conventional elas- 
tomers BSE [Hi. 

This novel behavior is a consequence of the semiffexible 
nature of filaments that make up the network. Such fila- 
ments are much stiffer than conventional polymer chains 
with persistence lengths Ip ranging from hundreds of 
nanometers to tens of microns [12|, • The persistence 
length is of a similar magnitude to other lengths in the 
system namely the total filament length L and the char- 
acteristic length between crosslinks or branch points Ic- 
As a consequence the filaments cannot be modeled as 
flexible chains {Ic ^ /p), nor as rigid rods {Ip ^ Ic)^ in- 
stead they are termed semiflexible because Ic ^ Ip. It is 
this interplay of three different length scales that distin- 
guishes the elastic behavior of biopolymer networks from 
rubbers where there is only one relevant length scale in 
the problem: the span between crosslinking points. 

There have been various theoretical ideas put forward 
to explain the elasticity of these semiffexible networks 




FIG. 1: A schematic illustration of the three relevant length 
scales in semiflexible networks. Because all three lengths are 
of a similar scale, some very interesting elastic behavior re- 
sults. 



[i, i, [Ti, [TBI, fl6|. In general these can be split into 
two classes. Models where the elasticity is explained in 
terms of the entropic stretching of individual filaments 
between crosslinks [8, 9, 14] and models where the elas- 
ticity is explained in terms of a transition from an en- 
thalpic bending-dominated regime to an entropic stretch- 
ing dominated regime [l5|, [ly]. One of the main dis- 
tinctions between the two approaches is that for mod- 
els considering only entropic stretching of filaments it is 
assumed that any deformation applied to the system is 
homogeneous down to the length scale of crosslinks, that 
is, the deformation field is affine [8]. For models relying 
on dynamic bending of filaments, the deformation must 
necessarily be non-afflne [HI, [13] • 

Recent experiments have claimed that the strain stiff- 
ening of such networks is satisfactorily explained in terms 
of the affine entropic stretching model [q, Q . In particu- 
lar Storm et al. [8] concluded that an affine model that 
accounts for the possible mechanical stretching of flla- 
ments is able to explain the elastic response of isotropic 
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fibrogen gels. This is in a curious contrast with Gardel et 
al. [2.], who reported that the elasticity of crosslinked F- 
actin filaments at various concentrations is explained in 
terms of entropic effects only, accounting for no mechan- 
ical stretching. In this paper we aim to investigate these 
ideas more closely by comparing an affine model predic- 
tion for the shear modulus of semiflexible networks to 
experimental data of both Storm et al. [8] (performed 
on the more flexible fibrogen networks) and also to the 
data of Gardel et al. [9] (performed on the more rigid 
F-actin networks). We also extend our model to the the 
case of uniaxial stretching, and compare it with the in- 
vivo results of Fernandez et al. [1^ who measured the 
stress-stiffening of entire cells under uniaxial strains. 

We develop a simple affine model that accounts 
for both entropic elasticity and the direct mechanical 
stretching of filaments in a similar way to Storm et al. Q . 
This is achieved by calculating the probability distribu- 
tion for the separations of crosslinks based on expressions 
obtained for semiflexible worm-like chain models [l9|, , 
which we extend to be valid at all extensions. Our model 
is therefore able to calculate the free energy density of the 
network explicitly and therefore we obtain any desired 
shear (or Young) modulus at a given imposed strain. We 
find that such a model does indeed explain the stiffening 
of more flexible biopolymer gels (such as fibrogen) rather 
well. This is in agreement with the model of Storm et al. 
0]. We also find that the same affine model is able to 
explain the stress stiffening of crosslinked networks of the 
stiffer filaments F-actin [9]. However, we find that if one 
properly accounts for the initial distribution of crosslink 
separations and orientations in the network, there is no 
longer a universal scaling of modulus with stress for the 
network as claimed in [9] . Moreover we find that the form 
of the modulus-stress relation becomes sensitive to the 
mechanical stretch modulus of the individual filaments. 
We also note a curious feature, that to achieve an agree- 
ment with the F-actin data we have to assume that the 
effective persistence length of the filaments in the net- 
work is much smaller than one usually finds in the lit- 
erature. In discussing this problem, we find that similar 
effect is known in various other systems, such as amyloid 
fibril networks, where many authors report a high persis- 
tence length on the basis of image analysis of curvature of 
filaments deposited on a substrate surface [21], while the 
few studies of bulk networks of these fibrils ^2,] give very 
much shorter persistence lengths. We propose that such 
a reduced persistence length might be a consequence of 
the filaments being quenched into a network. 



II. QUENCHED NETWORKS 

The quantities of interest in this paper are the equilib- 
rium elastic moduli of a crosslinked (or branched) semi- 
flexible polymer network. Any such equilibrium thermo- 
dynamic property is always expressible as a derivative of 
the free energy of the system, and so the task is really 




FIG. 2: Schematic illustration of a strand in a semiflexible 
crosslinked network. The contour length of filament between 
the two crosslinks is denoted by Ic while the vector connecting 
the two crosslinks is denoted by R. 

to be able to calculate the free energy of a crosslinked 
semiflexible network as a function of arbitrary imposed 
strain. The remainder of this section is devoted to this 
task. Let us initially consider a single section of a fila- 
ment that connects two crosslink nodes. We will refer 
to such filament sections as "strands". Each strand in 
the network is identified by the label a. The arc length 
of filament that makes up a strand is denoted by Ic and 
is assumed to be the same for all a (this is obviously a 
crude approximation, which however is expected to hold 
well in a large system with a relatively high crosslinking 
density). The displacement vector that connects the ends 
of a strand is denoted by Ra- For the remainder of this 
paper we will refer to the separation vector as the dimen- 
sionless vector r = R/lc These variables are shown in 
Fig. [2] 

The vectors will not all be the same, but will follow 
some distribution which is a characteristic of the fila- 
ments making up the strands. In general, we will label 
this normalized probability distribution by P{ra) - we 
delay finding an adequate physical expression for P{ra) 
until the next section. Given P{ra)^ the free energy of 
the network as a function of strain can be calculated in 
the following way: 

Firstly, we note that the crosslinks are quenched. That 
is, the time needed for a strand to explore different con- 
formations is much smaller than the time needed for 
the crosslinks to break and re-form. This separation 
of timescales means that one can perform equilibrium 
statistical mechanics over conformations while assuming 
that the topology of the network remains quenched in 
place - a classical procedure of quenched averaging em- 
ployed in theories of rubber elasticity [22|, 0, 0, [25[ . 

Secondly, we assume that the strands of the network 
do not interact apart from at crosslinks (i.e. they are not 
entangled on their own). This approximation is also well 
known from rubber elasticity and we expect it to hold 
for open- mesh semiflexible networks provided Ic/lp ^ 1- 
This is due to the fact that at such length scales the 
strands behave very rigidly and therefore have less op- 
portunity to explore the surrounding space, and hence 
we expect excluded volume effects to be negligible. Each 
strand a will have a free energy of: 

F(r„) = fcBTln-^ (1) 

If one now deforms the whole network by some strain ten- 
sor A, the separation vectors will individually trans- 
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form to a new separation which wih depend in some 
way on the imposed deformation: 



(2) 



where Aq;(A) is the strain tensor describing the transfor- 
mation of the ath strand. In general the strain tensor 
Xa win be a function both of the apphed strain A and of 
which strand we are considering, a - that is, the defor- 
mation does not in general have to be affine. 

As long as the topology of the network is conserved 
under the deformation, the free energy of the network 
can be expressed as the free energy of each strand in 
the new strained state F{Xara) summed over the same 
strand labels a: 



(3) 



We now make the assumption that the deformation 
field is affine. That is to say, that the strain tensors acting 
on each individual strand is the same as the bulk = 
A. We do not attempt to justify this assumption, but 
rather will examine the validity of it when comparing the 
results of this model to experimental results. With this 
affine assumption the free energy density of a quenched 
network undergoing affine deformation can be written as 
the integral: 



/(A) = nksT / d> P(r) In 
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P(Ar) 



(4) 



where n is the number density of strands. Provided one 
can obtain an expression for the probability of separation 
vector r occurring, one can calculate the free energy den- 
sity of the network using the above integral. The next 
section deals with the details of implementing this calcu- 
lation: in particular in finding an appropriate expression 
for P(r). 



III. THE MODEL 

This section will be concerned with developing an ex- 
pression for the probability distribution of semifiexible 
separation vectors P{r) that can be used in EqlH to 
calculate the free energy density of the network. The 
main difficulty is that although there are a number of ex- 
pressions that have been proposed for the form of P{r) 
HSSES. ah are calculated for inextensible chains 
and therefore have an essential zero at \r\ = 1. 

If one invokes the affine assumption however, some 
chains will inevitably be stretched beyond the limit \r\ = 
1. One therefore needs a probability distribution that 
accounts for the finite probability of the strand being 
in a mechanically stretched state - where the contour 
length of the strand has increased beyond Ic by means of 



stretching of bonds along the filament backbone. The 
work of Kierfeld et. al. [28] provides an analysis of 
single filaments that accounts for the possible mechan- 
ical stretching, however they succeed only in providing 
extension-force relations for such filaments in different 
force regimes, whereas we require the full probability dis- 
tribution. We achieve this in the following way: 

We again examine a single strand that is at its equilib- 
rium separation. When the separation of the strand ends 
r is at its equilibrium, we expect there to be no backbone 
stretching. The probability distribution is determined 
only by entropic effects of the semiffexible ff lament. In 
this regime P(r) is accurately described by distributions 
obtained in [l^, [2^, e.g.: 



P(r) (X 



(1 _ ^2)9/2 



Exp 



9a 



8(1 -r2) 



(5) 



where the parameter a = Ic/lp and is a measure of the 
ffexibility of the ff laments and r = \r\. As one stretches 
the strand from this equilibrium separation, the entropic 
returning force of the strand increases and it becomes 
increasingly difficult to stretch. The effective modulus of 
this entropic elasticity /ig is given by: 



p{r) 



(6) 



This is a function of extension r and diverges like (1 — r)^ 
in agreement with [29]. At some value of extension, 
which we will label ri, the entropic stretching modulus 
/ie(^) will exceed the mechanical stretching modulus of 
the strand jUrn (related to the Young modulus of the ffl- 
ament), /ie(ri, a) 

— I^m- This deffnes vi in terms of fiyn 
and a. We shall use a simple model for the cross-over 
from entropic stretching to mechanical stretching, which 
states that below ri the response of the chain is purely 
entropic, while above ri it is purely mechanical: 



Me(r) 



< r < ri 
ri < r < oo 



Entropic . . 

Mechanical ^ ^ 



Note we assume that the Young modulus of the fflament 
remains constant with strain. In reality it would decrease 
with strain until eventual rupture of the fflament, how- 
ever we do not deal with such terminal effects in this 
paper. From these expressions for the stretch modulus 
of an individual strand, we can obtain an expression for 
the free energy of the strand as a function of separation 
r by integrating twice: 



(3F(r) 



9a 



ln(l - H) + ^^^—^ 

2 * 

2 f-f^m^ 



< r < ri 
ri < r < oo 



(8) 



where A and B are constants that are determined by 
requiring that the free energy and its ffrst derivative with 
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FIG. 3: (a) P(r) plotted for values of a of 0.25 (solid line), 
0.5 (large dash), 1.0 (small dash), 5.0 (dots). Note that the 
last distribution (for more flexible chains) already resembles 
the Gaussian, with the peak at r = 0. The value of /Hm/kBT 
used here is 10^. (b) P(r) plotted for values of iim/kBT of 
5000 (solid hne) 1000 (large dash) and 500 (small dash). The 
value of a is 0.2 



respect to separation are continuous across the boundary 
at ri . Given this free energy we can obtain an expression 
for the probability distribution for the chain that is valid 
for all separations r using the usual definition: 



this will not in general occur at zero strain. Instead we 
must allow for the possibility that the network undergoes 
an initial spherically symmetric bulk volume change (an 
effect analogous to syneresis in ordinary gels). That is, 
we allow for the possibility that all strand separations un- 
dergo the transformation \r\ 6|r|, where 6 is a constant 
close to one. Only a spherically symmetric deformation 
is allowed since by symmetry it is the only spontaneous 
deformation possible in an isotropic network. We must 
therefore find the minimum of the integral 



47rr^dr P(r) \nP{hr) 



(11) 



with respect to h. We solve this numerically. As pointed 
out in [8], h is not an extra parameter in the problem, 
but is completely determined by a and firn- For net- 
works considered in this paper, all values of h are less 
than one; the network therefore shrinks initially after 
crosslinking, which is a well-known effect of reducing 
the conformational freedom of filaments on establishing 
quenched crosslinks. Once the minimum of the free en- 
ergy has been found, we calculate the free energy density 
as a function of strain using EqHl We do this for two 
deformation geometries: shear and uniaxial extension. 
These are implemented as follows: 



A. Simple shear deformation 

The experiments [l, [9[ impose a simple volume- 
conserving shear deformation on the networks. Trans- 
formation of the separation vector components (x, z) 
under simple shear depends on a single parameter 7: 



P(r) = 

where jS = l/ksT and 



-0[F(r)-F] 



(9) 



(10) 



We remind the reader that the separation is a vector, and 
therefore integration must be performed over all orienta- 
tions as well as the length of r. Plots of P(r) for various 
values of the parameters a and fim are shown in Fig. [31(a) 
and (b). 

The advantage of our approach is now clear; we have a 
relatively simple expression for the probability distribu- 
tion of strand separation vectors that is valid for all sepa- 
rations, from which we can calculate any thermodynamic 
property related to a single strand. With the probability 
distribution defined, we can proceed to calculate the free 
energy density of the network as a function of applied 
strain via EqlH This is implemented as follows: 

Firstly we must find the initial minimum of the free 
energy as this is the equilibrium point where the net- 
work starts before being strained. As observed in 





(12) 



Transforming to spherical polar coordinates we can 
write the norm of the transformed separation in terms 
of the original separation r, the angles parameterizing the 
orientation of a strand (^, (j)) (defined in the usual way 
with respect to the polar axis z) and the applied shear 7: 



r' = hr^J 1 + 27 sin d cos d cos + 7^ cos^ d (13) 
The integral to be evaluated is therefore: 



/(7) 



-nfcsT y" P(r) In P(r')r2 sin Mr d6' d(/) 



(14) 



This integral is an expression for the free energy of a 
semiflexible network that has undergone affine deforma- 
tion at any shear strain 7. The integral is performed 
easily in Mathematica. There is sometimes a confusion 
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in the hterature about which quantity to cah the "elas- 
tic modulus" when the deformation is not infinitesimal. 
Two different definitions are possible: the shear modulus 
G and the differential shear modulus K of the network 
that has been subject to a simple shear of magnitude 7 
in the xz plane are easily obtained from the free energy 
density by using the following relations [s^]: 



G 



2dl 

7 dj 



K 



dcr^r, 



^ 2 — 

97 ^7^ 



B. Uniaxial extension 

There are a number of experiments that impose uni- 
axial compression or extension on the cytoskeleton, in 
vitro ^0] or in actual living cells [18]. To compare this 
afiine model to such experiments we must calculate the 
modulus in this geometry. Assuming this deformation is 
volume conserving, the affine transformation of the sep- 
aration vector is: 



1/VA 
1/VX 
OA 




(16) 



As before, expressing this in spherical polar coordi- 
nates, the norm of the separation vector transforms as: 



br\ 



' sm 



A 



A^ cos^ 



(17) 



Following the same arguments as in the case of sim- 
ple shear, the free energy density of a network that has 
been subject to affine uniaxial volume conserving strain 
is therefore: 



/(A) 



-2TTnkBT J P{r) hiP{r')r^ sin Mr (16*, 



(18) 



where we have picked up a factor of 27r from the blank 
integration over 4>- Given this elastic free energy density, 
the Young modulus Y and differential Young modulus E 
of the network are given by: 



Y 



idl 

XdX 



d(Tz 



(19) 



is how the linear modulus of the network (at infinites- 
imal deformations) scales with the parameter a = Ic/lp, 
and how the non-linear modulus scales with applied pre- 
stress. For the affine model considered in this work we 
find that the linear shear modulus Go and the linear ex- 
tension modulus Yq have a scaling with a of the form: 



(15) ^0 ^ LlnksT 



Yq ^ 3.3nkBT 



(20) 



Naturally, both linear moduli are independent of /i^ at 
infinitesimal strains, and they are related by a factor of 3 
as indeed required in an incompressible linear elasticity 
[30] . The scaling of is how the "modulus" of sin- 
gle filament scales for large strains. The contribution of 
the form (4 — a)^ comes from the quenched averaging. 
Once we quench the network with cross-links, we allow a 
bulk spherically symmetric transformation of all separa- 
tion vectors, to minimize the quenched free energy. This 
contribution accounts for the (4 — a)^ dependence. Note 
that for small a = Ic/lp (very rigid filaments), the scaling 
is approximately inverse square. 

It has been claimed that for models considering the 
simple affine entropic stretching of constituent strands, 
the scaling of the modulus K with stress a should follow 
a simple power law K{a) ~ a^/^ We find that this 
reasoning is, however, only valid if one considers a single 
chain. If instead one properly treats the orientational av- 
eraging over the distribution of separation vectors, such 
universal scaling is lost. We illustrate this effect in Figs. 
[H and [5] by plotting the stiffness of a single filament as 
a function of tension, for various values of the stretch 
modulus jiim at fixed a. In Fig. [H we plot the "modu- 
lus" of a single chain dr/dx^ where r is the tension on 
the filament. As expected, at high r the filaments show 
a power law scaling ~ r^/^: it is indeed true that for a 
worm-like chain model of a single semi- flexible filament, 
the modulus scales like the tension to the power 1.5 for 
large extensions, exactly as argued in [9] and [31^. How- 
ever, if we have a crosslinked network with a collection of 
filaments with a corresponding distribution of separation 
vectors P(r), this result changes. In Fig. [5] we plot the 
"modulus" d{r)p/dX as a function of "stress" (r)p for 
different value of /i^n, where {..)p refers to averages over 
the probability distribution of separation vectors. We 
find that no such universal scaling exists - the functional 
dependence is no longer a simple power law. In particu- 
lar, the initial non-linear exponent can be made greater 
than 3/2 provided /i^n is large enough. 



IV. RESULTS 



C. Scaling of moduli 

For both shear and uniaxial-extension geometries, we 
can examine how the stiffness of semiflexible networks 
scales with model parameters. Of particular interest 



We now compare the results of the model presented 
above with a series of experiments. Firstly we examine 
the data of Storm et al. [8] , who measure the simple shear 
modulus G of branched fibrogen networks as a function 
of applied strain at various fibrogen concentrations. The 
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FIG. 4: Log-log graph of the modulus of a single filament as 
a function of applied tension, where there is no initial dis- 
tribution of separation vectors. The curves correspond to 
= 10^, 10^ and 10^ (top to bottom), for a fixed value of 
a = 1. The constant gradient at high tensions corresponds to 
a power law scaling of dr/dX ^ r^^^) 
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FIG. 5: Log- log graph of the average modulus as a function 
of the average tension, where averages are performed over 
the distribution of separation vectors P(r). Curves are for 



10^ 



10^, 10^, 10^ (top to bottom). Accounting for 



this distribution destroys the simple power law scaling of mod- 
ulus with tension, even in a 1 — d network. The dashed line 
shows a power law of 3/2. 



data of [8] is taken over a wide range of strains (up to 
~ 100%) and allows the model to be tested well into the 
regime of non-linear elasticity, up to the point of network 
rupture. It does not, however, provide a good opportu- 
nity to examine the behavior of such networks at low 
strains (< 5%) because the first data point already oc- 
curs at 5% strain. 

Secondly we compare the predictions of our affine 
model to the data of Gardel et al. [9], who examine the 
differential modulus K of crosslinked actin networks, as a 
function of applied prestress at various concentrations of 
actin monomer, while keeping the ratio of monomer con- 
centration to crosslinker concentration constant. This 
enables us to examine more closely how well an affine 
model compares to experimental results in stiffer net- 
works at smaller strains. 

We also compare our model to the experimental re- 
sults of Fernandez et al.[18] who measure the uniaxial 
stiffness of entire living cells as a function of pre-stress. 
They find that the uniaxial stiffness exhibits a universal 
scaling with applied pre-stress, which we find is consis- 
tent with cells, whose stiffness is determined by an actin 




FIG. 6: Schematic illustration of many strands of length Ic 
along a filament of length L. The approximate number of 
such strands on a single filament is m ~ L/lc. 



cytoskeleton. 



A. Comparing to data of Storm et. aL 

Before we compare our model to experiment we must 
choose how we are to fit to data. The model we formu- 
lated in EqHlhas three parameters: n (the number den- 
sity of strands in the network), a (the ratio of the length 
of the strand to the filament persistence length Ic/lp) and 
lira (the mechanical stretching modulus of the filaments). 
However, the number density of strands n and the value 
of a are not independent. They are related in the fol- 
lowing way: Let us say that we have a filament of total 
length L. Crosslinks or branching points are spaced by an 
average length Ic along the filament backbone. Assuming 
Ic <^ the number of strands m on a single filament is: 
m ^ L/lc. This does not account for dangling ends, but 
in a densely cross-lined network this becomes irrelevant 
as m and m — 2 would differ by a small fraction. These 
lengths are shown in Fig. [6l 

The total number of such filaments in our system must 
approximately be = Ltot/L where Ltot is the total 
length of all filaments in the system. This total length of 
all filaments can be written as the total mass of filament 
monomer, M, divided by the mass per unit length of fil- 
ament, p. So the number density of filaments is Cm/Lp^ 
where Cm is the mass density of monomers, Cm = M/V . 
Therefore the total number density of strands in the sys- 
tem is: 



Lp 



-m ■■ 



Icp 



alpp 



(21) 



Since Cm, lp and p are independently measured exper- 
imentally, the value of n is determined directly by the 
value of a. In examining experimental data, we will fit a 
and n independently, and then use the above relation as 
a check on consistency of the fits. We propose that if the 
fitted parameters a and n are consistent, then the value 
of nalpp/cm should be close to unity. 

The parameter is defined above as the energy scale 
associated with a fractional increase in filament length. 
We can obtain an estimate of the value of p^yi 

for a strand 

by realizing that the energy of a filament with a Young 
modulus Ef and cross-sectional area A that has been 
stretched from a length Ic to / is: 
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FIG. 7: (a) Fit of the model to the data of Storm et al. 
0]. The data sets are obtained from shear measurements 
on branched fibrogen networks with monomer mass concen- 
trations of 4.5 (♦), 2.0 (■), 1.0 (•) and 0.5 (□) mg/ml. The 
model fits very well to the data at all concentrations for phys- 
ically reasonable and consistent fit parameters fTable ITV A|) . 
although fails to capture the concave nature of the curves at 
high strains, thought to be due to filament rupture, (b) The 
same fit as in shown on a Linear- linear scale. 



8 = ]^Efl,A{l/l, - If (22) 
Comparing to the original expression for /i^ we see: 

Mm EfAlc EfAlpa (23) 

An estimate for the mechanical stretching modulus can 
therefore be obtained, provided one knows the Young 
modulus of the filament and the approximate radius of 
the filament. In comparing our model to the data of 
Storm et al. we fix the value of /i^n to be bOOOksTa. 
We do this because the Young modulus of the filaments 
should be the same at all concentrations. This is con- 
sistent with the measured diameter of lOnm for fibrogen 
filaments and a Young modulus of ~ IMPa, which has 
been measured experimentally for fibrin [32]. In that 
case, only one truly independent fitting parameter re- 
mains (the stiffness ratio a) and thus the test is rather 
stringent. The resulting fits are shown in Fig. [71(a) and 
(b), and the fitted values of a are given in Table. (|lVAj). 

Note that the check of consistency for the fitted values 
all fall close to the estimate value of 1. We believe this 
shows that the affine assumption is valid for the above 
kinds of networks. 



Cm{mg/ml) 4.5 2.0 1.0 0.5 

fim/kBTa 5000 5000 5000 5000 

a 0.60 0.75 0.95 1.40 

b 0.994 0.992 0.987 0.968 

kBTn{J/m^) 9.7 2.2 0.69 0.32 

nalpp/cm 1.89 1.21 0.96 1.31 

TABLE I: Parameter values for the fits to data of Storm et al. 
[1] . These fit values have been obtained using the persistence 
length Ip of 0.5/im (reported in ^]), and a mass per unit length 
for fibrogen p of 2.4 x lO-^'^kg/m [|. 



B. Comparing to the data of Gardel et aL 

To compare the affine model proposed here to the data 
of Gardel et al. [9], we use a slightly altered method to 
that described above. We do this for the following rea- 
son: Gardel et al. perform measurements of differential 
shear modulus K on scruin-crosslinked F-actin filaments. 
They measure both the number concentration of actin 
monomer in the system ua as well as the number con- 
centration of scruin crosslinker in the system n^. They 
then measure the differential shear modulus as a function 
of pre-stress for four values of ua while keeping the ratio 
u = Us/riA = 0.03 constant. Recalling Eq|2T] we have 
a relationship between the number density of strands 
in our system n and the contour length of the strands 
Ic'. n = Cm/lcpi where Cm is the mass concentration of 
monomer and p is the mass per unit length for F-actin 
filaments. If such strands are formed at crosslinks, and 
each crosslink has an average coordination of z (we shall 
assume 2: = 4 here), then the number density of strands 
n is related to the number density of cross-links Uxi via 
n = nxiz/2. Therefore knowing Uxi from experiment, we 
have an estimate for Ic via: 

/c = ^ = ^ 0.04/.m, (24) 
pzuxi 2pu 

where in the second term we have expressed Ic in terms 
of the mass of a single actin monomer rriA (taken to be 
42kD (33I, Isj] ) and the ratio of crosslinker concentration 
to monomer concentration v = Us/ua which was kept at 
a constant value of 0.03 in these experiments. In the final 
term we have evaluated the contour length of a strand for 
the experimental values of [9] and p = 2.6 x 10~^^kg/m 
[3^ . Note that this predicts that Ic is only dependent 
on the ratio of crosslinker to actin monomer v = Us/tia 
and so is constant for the different data sets of Gardel 
et al. This is what we should expect: if strands are only 
formed at crosslinks as we suppose, then doubling both 
the number of crosslinks and the number of monomers 
should not result in a change of Ic. This is in contrast 
with the argument offered in |9|] , where the authors pro- 
pose that Ic is proportional to the entanglement length 
(in turn, correctly estimated in [35]) at a fixed u. 

Having an estimate for Ic we now estimate the value 
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Cm{mg/ml) 1.2 0.90 0.50 0.35 

^im/kBT 6 X 10^ 6 X 10^ 6 X 10^ 6 x 10^ 

a 0.005 0.005 0.005 0.005 

kBTn{J/m^) 4.4 3.2 1.8 1.2 

TABLE II: A table of the model parameters we would expect 
if we take the persistence length of filaments to be Ip — lO/xm. 
Such a value for Ip would result in very stiff filaments. From 
Fig. [8] (a) we see that such model parameters produce K{a) 
curves that drastically fail to agree with experimental data. 



of the mechanical stretching modulus firn from the ex- 
pression jirn = EfAlc, where Ef is the Young modulus 
of F-actin filaments and A their effective cross-sectional 
area. We can estimate the value of fim for such filaments 
from the value for EfA = 5 x 10~^N measured for F-actin 
in [36]. This yields jUm = 6 x 10^ ksT. Since the value of 
Ip is constant for the data sets, the value of jUm niust also 
be constant. If we now take a value for the persistence 
length of actin of 10/im which is most frequently quoted 
in the literature [12, 37, 38], we then have estimates for 
the three parameters of our model (n, /i^n, taken di- 
rectly from experimentally measured quantities. These 
are shown in Table [Til 

The predicted curves of differential modulus against 
applied pre-stress for the parameter values in Table HIl are 
shown in Fig. [8] (a) along with the plots of the original 
data. It is clear such parameter values give very poor fits 
to data predicting an unstrained modulus K{0) that is 
out three orders of magnitude. 

We would like to point out that in order to get a good 
fit to experimental data it is not simply a case of rescaling 
the number density of filament strands in the network 
n. Indeed doing so gives curves that fail to predict the 
onset of stress stiffening by over an order of magnitude. 
Remarkably however, we do find that the affine model 
can provide very good fits to the data of [9] as shown in 
Fig. [8] (b), for the different set of fitting parameters in 
Table. Illli There are a few interesting points to note. 
We have fixed the value of both /i^ and Ic for such fits 
as we have estimates for both, based on experimentally 
measurable quantities. We then perform a two-parameter 
fit by varying Ip and n. Good fits can only be achieved 
with this afiine model by using a value for the ratio a = 
Ic/lp, that is two orders of magnitude larger than would 
be expected from the simple arguments given above - 
that is, a value of Ip = 0.022/im. As before, we check 
the consistency of the values for n by calculating nlcp/cm 
which should be close to unity - and indeed we find values 
that are reasonably consistent. Note we do observe a 
scaling resembling K{a) ~ a^/^, in agreement with the 
statement in [9]. However, we find that such scaling is 
not universal for such networks but depends strongly on 
the values of a = Ic/lp and /irn- It is in fact because of 
the particular fitted values of a and /irn, that we observe 
a scaling like K{a) ~ a^/^. 

In short, we find that assuming the persistence length 
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FIG. 8: (a) K{a) predicted by the affine model for the param- 
eters in Table |ll] compared to the experimental data. If one 
assumes that Ip = lO/xm, then the affine model produces very 
poor agreement with the data. Apart from the disagreement 
in values, the model does not show the scaling of K{a) ^ cr^^^ 
observed, (b) Two-parameter fits of the model to the data 
of Gardel et al. ^ on log-log axes. The data shown corre- 
sponds to differential shear modulus for networks with F-actin 
monomer concentrations ca of 29.4 (□), 21.4 (•), 11.9 (■) and 
8.33/iM (♦) and constant value u — 0.03. The fitted value 
for Ip = 0.022/im, assuming the value of Ic remains fixed at 
0.04/im. The solid straight line shows a scaling of a^^^ . 



of F-actin filaments in these networks is ~ 10/im results 
in the affine model failing severely and predicting that 
the networks would be far stiffer than is actually mea- 
sured. More importantly, there is no easy way out of this 
disagreement because experimental data have several key 
features that are all linked. However, we find the affine 
model of filament network can be reconciled with experi- 
ment, and produce very good agreement provided we as- 
sume that the F-actin filaments behave as if they are far 
more flexible than is usually quoted - with a persistence 
length of ~ 0.02/im. Interestingly such behavior has ac- 
tually been reported experimentally in networks of actin 
fllaments, who report values of the persistence length as 
low as 0.5/im [33|, ll^ which, although still much larger 
than the value needed by the affine model, shows that the 
persistence length of filaments embedded in networks can 
be far less than what is most frequently reported for sin- 
gle actin filaments. It is, of course, also possible that 
the in- vitro simple shear experimental geometry at large 
pre-stress underestimates the actual modulus, so that Ip 
closer to 0.5/im would actually be acceptable. 
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Cm{mg/ml) 


1.2 


0.90 


0.50 


0.35 


l^m/kBT 


6 X 10"^ 


6 X 10^ 


6 X 10^ 


6 X 10 


a 


2.0 


2.0 


2.0 


2.0 


b 


0.92 


0.92 


0.92 


0.92 


kBTn 


3.75 


1.25 


0.52 


0.29 


TllcX/ Cm 


0.89 


0.40 


0.30 


0.24 



TABLE III: Model parameters for the fit to the data of Gardel 
et al. shown in Fig. [8](b). It is a two-parameter fit for the 
values of n and a. Fitted values are obtained assuming the 
contour length of Ic = 0.044/im and a mass per unit length for 
F-actin p = 2.6 x 10~^^kg/m. The fitted persistence length is 
therefore 0.022/im. Provided we take the persistence length 
to be so small, we achieve good fits to data with reasonably 
consistent values for n. 



E/Eo 




FIG. 9: Fit of the model to the collapsed data of Fernandez 
et. al. [18]. A very good fit to data can be achieved by the 
parameters a = 0.5 and ppm = 5 x 10^. These values are 
consistent with an act in network. 



C. Comparing to the results of Fernandez et. al. 

To highlight the fact that this model might well be of 
biological significance, we compare its predictions with 
the data of Fernandez et al. [18] who measure the uni- 
axial stiffness E of entire cells as a function of applied 
pre-stress. They have found that the stiffness of a cell 
follows a simple master equation: 

^(r) = ( ^ = ^0 < (25) 

which is to say, a collapse of the modulus stress curves 
is achieved by rescaling the modulus axis by £^0, and 
rescaling the stress axis by dc. The model presented here 
provides a possible explanation for this behavior. The 
non-linear behavior is a consequence of chains being ini- 
tially entropically and then subsequently mechanically 
stretched. The scaling behavior of E{a) is a function 
of our key parameters a and /i^, that is, it depends 
on the transition from entropic stretching to mechani- 
cal stretching. We fit this data by for a value of a = 0.5 
and /j^m = 5 x 10^/cbT, with the result shown in Fig. 
[9l The value of the fitted parameter = 5 x 10^ is 
consistent with filaments with Young modulus Ef ~ 10^ 
Pa, diameters ofd ^ 5nm and lengths Ic ~ 0.5/im. These 
values are all consistent with an actin network. This sug- 
gests that the stiffness of such cells could be determined 
by the actin cytoskeleton network, with strand lengths 
of ^ 0.5/im and persistence lengths of ~ Ifim (note that 
this is also much smaller than the 'accepted' values of 
Ip ~ 10/im). 



V. CONCLUSIONS 

In this work we have shown that the equilibrium shear 
modulus of branched fibrogen networks and crosslinked 
F-actin networks as a function of applied strain is satis- 
factorily explained by an affine model that accounts for 
mechanical as well as entropic stretching. In the case 



of fibrogen networks we find the affine model fits the 
data remarkably well with fit parameters which are both 
consistent with each other and with experimental mea- 
surements on individual fibrogen fibers. In the case of 
F-actin networks we find that the model predicts the ob- 
served stress stiffening behavior accurately provided we 
allow the effective persistence length of the filaments in 
the network to be more flexible than is expected. That 
is to say, an affine model agrees with the data well and 
has reasonably consistent fit parameters provided the ef- 
fective persistence length of filaments in the network is 
taken to be Ip ~ 0.02/im. This is far lower than the ob- 
served persistence length of individual F-actin filaments 
reported in individual measurements [11,, ^37,], but a re- 
duced persistence length {Ip ~ 0.5 /im) has been rep orted 
experimentally in some F-actin networks [SS*, 39], and 
also arises from the fitting of the in- vivo data of Fernan- 
dez et al. [iH. 

This suggests to us that an afiine model does rather 
well at explaining the elasticity of branched or crosslinked 
semi-flexible filaments because in such networks the fila- 
ments behave far more flexibly than expected on the basis 
of measurements of filament curvature, when they are ad- 
sorbed on a substrate (note that the few measurements 
which measure Ip with quasi-elastic light scattering of fil- 
aments in networks often report much shorter Ip [33, 39|). 
Because the length of filament strands, /c, is actually of a 
similar order or even greater than the effective persistence 
length of the filaments, treating the strands as indepen- 
dent is a reasonable approximation - as it is for rubber 
- where one does indeed expect the affine model to work 
well. For more rigid networks that are composed of fil- 
aments where the effective persistence length is smaller 
than the contour length of strands, Ip Ic, we would 
expect the affine assumption not to be valid, and mod- 
els incorporating non-affine deformations may well need 
to be considered [HI, [l7|, SO]. However, it appears that 
in- vitro actin networks extensively studied in the recent 
years are not in this regime. 

We also find that properly accounting for the distri- 
bution of cross-link separation vectors has a profound 
effect on the scaling of the modulus of the network with 
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stress. Averaging over the distribution of initial separa- 
tion vectors means we do not see a universal scaling of 
the modulus with stress of the form K{cr) ~ a^/^ as pre- 
dicted by [9], but rather a functional dependence that is 
not a power law, and whose form depends on the values of 
a = Ic/lp and fim for filaments in the network. The fact 
that such a scaling is non-universal, but depends on chain 
parameters, can be used to explain why different scal- 
ing regimes are observed in different networks. Although 
Gardel et. al. [9] observe scaling like K{cr) ~ a^/^, ex- 
periments on dendritic actin [l^ have observed a scaling 
like E{a) ~ <j^"^, while experiments on entire cells [18[ 
have observed a scaling like E{(j) ~ a. We have shown 
that an affine model is able to capture the different scal- 
ing behaviors observed by accounting for the mechanical 
stretching of filaments. This could be of biological sig- 



nificance, as the non-linear behavior of the modulus will 
in general depend on the parameters of the constituent 
filaments. This suggests that one can tailor the form of 
the non-linear behavior of such networks by tuning the 
filament parameters a = Ic/lp and /i^ = EfAlc. That is, 
the form of the non-linear modulus as a function of stress 
K{a) can be drastically altered by changing Ic- 
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